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An analytical approach is developed to the problem of computation of monotone Riemannian 
metrics (e.g. Bogoliubov-Kubo-Mori, Bures, Chernoff, etc.) on the set of quantum states. The 
obtained expressions originate from the Morozova, Cencov and Petz correspondence of monotone 
metrics to operator monotone functions. The used mathematical technique provides analytical 
expansions in terms of the thermodynamic mean values of iterated (nested) commutators of a model 
Hamiltonian T with the operator S involved through the control parameter h. Due to the sum rules 
for the frequency moments of the dynamic structure factor new presentations for the monotone 
Riemannian metrics are obtained. Particularly, relations between any monotone Riemannian metric 
and the usual thermodynamic susceptibility or the variance of the operator S are discussed. If the 
symmetry properties of the Hamiltonian are given in terms of generators of some Lie algebra, the 
obtained expansions may be evaluated in a closed form. These issues are tested on a class of model 
systems studied in condensed matter physics. 
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I. INTRODUCTION 


1. The notion of distance has found 
[0, geometrical description of phase 


In recent decade different concepts emerging from quantum-information geometry have been intensively incor¬ 
porated in condensed matter physics |5l-ll7l|. The underlying idea may be briefly reviewed as follows. The properties 
of macroscopic phases of matter should be encoded in the structure of the quantum state. It is defined by a density 
matrix p depending via the Hamiltonian 'H on a set of parameters, e.g. coupling constants A and external ’’fields” 
h. If 'H smoothly depends on the coupling constants and external fields one has a map between (A, h, the inverse 
temperature /3) and the set of density matrices: p = p{P, A, h) := [2'(/3, A, h)]~^ exp[—/3H(A, h)], where Z(/3, A, h) is the 
partition function of the system. Further, the set of density matrices p should be endowed with a metric structure and 
thus explored as a Riemannian manifold. The realization of the set of density matrices p as a Riemannian manifold 
implies to define a relevant distance between two mixed quantum states [l|, 
applications in different fields. For example: thermodynamic of small systems 
transitions and quantum criticality 
testing and discrimination of states 


quantum estimation of Hamiltonian parameters [18l Il9l| and hypothesis 
, are only a part of them. If we consider the distance between two states 
obtained by an infinitesimal change in the values of the parameters that specifie the quantum state, we come to the 
notion of a metric tensor, i.e. the set of coefficients of the linear element ds^ when written as a quadratic form in the 
differentials of these parameters. In the present study the infinitesimal variation of a single parameter, e.g. the field 
h is considered, and the notations of a linear element and metric will be used interchangeably. 

For our purposes we shall recall some metrics which are interesting in the context of their physical applications. 

The metric after the name of Bogoliubov, Kubo and Mori [23 - [^ . has a clear physical meaning since it 

describes the isothermal susceptibility of the system. Along the years the Bogoliubov-Kubo-Mori (BKM) metric (under 
different names: Kubo-Mori scalar product, Bogoliubov inner product, Duhamel two-point function or canonical 
correlation) has been intensively studied in the community of both physicists and mathematicians , see e.g. 

IsBjl and references therein. 

Notably, the Bures metric depending on the fidelity between two density matrices has obtained increasing 
popularity in the information approach to physics, see e.g. [1, [13, and references therein. The Bures metric 

appears under different names: quantum Fisher information (apart from a numerical factor), SLD metric, fidelity 
susceptibility, etc [3.13ilsL[l3|. 

Much as for quantification of the asymptotic behavior of the error in the quantum state discrimination problem, 
one may define the quantum Chernoff metric dgQ, which is expressed in terms of the non-logarithmic variety of the 
quantum Chernoff bound|2lJ-|23| . The quantum Chernoff metric also enjoys increasing popularity in the information 
approach in physics It appears under several different names, among which are (up to a factor of one forth) 

the Hellinger metric [J or the Wigner-Yanase metric [2ll - l^ . 

Another Riemannian metric in the physical setting, (the origin of the subscript MC will become clear in 

the subsequent consideration), is related to the quadratic fluctuations of a quantum observable (the variance of this 
observable). The metric was obtained as an approximated form of the dispersion of the Umegaki’s relative 
entropy or as the Hessian of the quasi-entropy considered in [^ . 

Here we stress that all the mentioned above metrics d^^^, d\, dq^ and d\jQ belong to the large class of the 

monotone (or contractive) metrics. In the geometrical approach to statistics proposed by Morozova and Cencovp^ 
the monotone metrics can be introduced and studied from a unified point of view due to the work of Petz [d^ • 

In view of Morozova, Cencov and Petz studies it is established that there exists one to one correspondence between 
monotone metrics dj and Ldvner operator monotone functions ” /” (on the subject of the operator monotone functions 
see, e.g. the monograph ^). Each one of the metrics may be written as a sum of two contributions; the classical 
Fisher-Rao term and a quantum term that depends on the definition of the inner product that induces the monotone 
metric on the space of quantum states. The last is far to provide uniqueness. The consequences of this fact are 
subject of considerable interest, see, e.g. along with a number of references therein. The relation between 

operator monotone functions and monotone metrics allows to explorer in deep the relations between the different 
metrics in the context of quantum statistical mechanics and condensed matter physics. 

The aim of the present study is to show that the isothermal susceptibility (respectively d^^^) and the variance 
of the operator involved through the control parameter (respectively d^^) may serve as reference metrics for the 
whole class of monotone metrics. In other words, each one of the monotone metrics can be divided into two parts: 
the d^p^^ (or d\jQ) and a constituent in the form of an infinite series of thermodynamic mean values of products 
of iterated commutators which are related with the noncommutativity of the problem. It is shown that this suggests 
a neat and deep interplay with the linear response theory through the relation with the moments of the Dynamical 
Structure factor (DSF). If the Hamiltonian is a linear form of the generators of a Lie algebra the final result may be 
presented in a closed form. 

The paper is organized as follows: In Sec.H, using the fact that all monotone Riemannian metrics are characterized 
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by means of operator monotone functions we derive a generic spectral representation of the monotone Riemannian 
metrics (see Ea. dT^ '). We show that this formula provides a sufficiently simple way to derive different well known 
metrics, e.g. Bogoliubov-Kubo-Mory, Bures, Chernoff and Morozova-Cencov metrics. Using this spectral represen¬ 
tation, in Sec. Ill we obtain new series expansions for the monotone Riemmenian metrics [see Eqs. m and (EH)]- 
In Sec. IV, the proposed series expansions are checked against explicit expressions in a concrete model. A summary 
and discussion are given in Sec. V. The Appendix is a compendium of known and not so well known facts about 
the Bogoliubov-Duhamel inner product, the functionals introduced in [sslj and their relations to the moments of the 
dynamical structure factor. 


II. MONOTONE RIEMANNIAN METRICS: GENERIC FORMULA 

Following the work of Petz [dj (see also Si) let us recall that a Riemannian metric on the manifold D of the 
density matrices p can be written in the form 

gpiX,Y) = {X,mf{Lp,Rp)-\Y))HS, (1) 

where X, Y belong to the tangent space TpV at p of the manifold 2?, Lp{A) := pA and Rp{A) := Ap stand for the left 
and right multiplication by p for A belonging to TpV, 

mf{A, B) ■= (2) 

is the Kubo - Ando operator means and {AB)hs '■= Tr{A*B) is the Hilbert-Schmidt inner product. Hereafter the 
subscript / stands for a monotone metric which depends on the operator monotone function f{x) : = f{x)/x 

with /(I) = 1. A Riemannian metric Pp defines a square distance cPj between two infinitesimally close states p and 
p + dp which in the basis independent form is given by 

d) = gp{dp,dp). (3) 

Note that since the operators Lp and Rp commute, the following relation holds mf{Lp,Rp)~^ = Cf{Lp,Rp), where 
the symmetric function Cf{x, y) = Cf{y,x), is called the Morozova-Cencov function. Recall that the Morozova-Cencov 
function has the following presentation 


Cf{x,y) = 


xfiy/x)' 


The function Cf{pm,Pn) obeys the equation Cf{tx,ty) = t ^Cf{x,y) for any t G R. 

From now for concreteness we shall consider one-parameter family of density matrices 

p{h) = [ZAr(/i)]“^ exp[-/3?2(/i)], 

defined on the family of iV-particles Hamiltonians of the form 


(4) 


(5) 


nih) = T-hS, 


( 6 ) 


where the Hermitian operators T and S do not commute in the general case. Here, /i is a real (control) parameter, 
ZN{h) — Trexp[—/372(h)] is the corresponding partition function and /3 = {KbT)~^ is the inverse temperature. We 
assume that the Hermitian operator T has a complete orthonormal set of eigenvectors jm) with a non-degenerate 
spectrum {Pm}; T\m) = Tm\m), where m = 1, 2,... In this basis the zero-field density matrix p := p{0) is diagonal: 

{m\p{0)\n) = pmSm,n, Pm := e“^'^™/^Ar(0), m,n=l,2,.... (7) 

Now, let us consider two nearby states pi = p and p 2 = p + dp produced by infinitesimal changing of the control 
parameter h. We choose p = p(0) and consider the matrix elements of dp, e.g. {m\dp\n), in the basis where p is 
diagonal. Thus, for p = p(0) with matrix elements ([7]) through ([3]) any monotone Riemannian metric on the set of 
quantum states (up to a proportionality constant) is presented explicitly as (see, e.g. Mi 13): 



771 


dpL 


\ XI Cf{pm,Pn)\{m\dp\n)\'^. 
m,n,m7^n 


Pm 


( 8 ) 
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It is however possible to give a more convenient (for calculations) presentation of by adopting the statistical 
mechanics viewpoint which will be presented below. 

Following @ (see also for the first term in Eq.([8]) we have 


1 

4 


E 


dpL 

Pm 



(9) 


where 5'^ is the diagonal part of the operator S and 


:= [Z(T)]-iTr{exp[-/3T]-..} 


( 10 ) 


denotes the thermodynamic mean value. A next step consists in using the following relations (see, e. 


g.S[Tl): 


{m\dp\n) = {m\dn){pn - Ptti), m ^ n, 


( 11 ) 


and 


mldn) = (’yWI") , mi.) = ^ = -S, 


Tn-T^ ’ dh 

By means of Eqs. ®, (HH) and m, Eq.® takes the more convenient form 




{{ 5 S'^f)T+ Y. ^fiPmiPn) (inpEln™™) 

m,n,m^n 


\{m\S\n)\‘^ 


( 12 ) 


(13) 


This expression is a basic element in our further calculations. 

There are so many monotone metrics as /’s are. Ones that are mostly proposed in the literature operator 
monotone functions are as follows [d^: 


fnarix) = 


2x 


, fsKAlix) = 


X — 1 
Inx 


,fwY{x) = 


(1 + y/xf 


Jb{x) = 


1 


, fMc{x) = 


X — 1 
Incc 


1 + a; 


(14) 


The function fnarix) is the minimal, while fsix) is the maximal operator monotone functions on [0, +oo). The former 
defines a metric known as RED metric while the last one gives rise to the SLD metric (named also Bures metric or 
fidelity susceptibility). The function fsKMix) leads to the Bogoliubov-Kubo-Mori metric. The function fwY{x) is 
associated to the Wigner-Yanase metric (or an equal ground the name quantum Chernoff metric is reasonable ). The 
function funix) was first conjectured in the paper |43j . which explains the subscript MC. Its matrix monotonicity 
was proved in |44| . 

It is instructive to consider the generic formula Eq. m in the concrete cases of the operator monotone functions 

(HI- 


A. Bogoliubov - Kubo - Mori metric 


The derivation of the BKM metric and its physical meaning as susceptibility ware discussed earlier in [T3, . 

In order to introduce the needed definitions we shall give our slightly different derivation based on the notion of 
the Bogoliubov-Duhamel (or Kubo-Mori) inner product (see the Appendix). The corresponding Morozova-Cencov 
function (everywhere below we shall use the subscript of the function / instead of the very function) is 

f\npr,-\npm\ 

CBKM\Pm,Pn) = - ■ (15) 

\ Pn Pm ) 

Erom the other hand, let us consider the spectral representation of the Bogoliubov-Duhamel inner product [2^ [H, 
mi) nil (see also the Appendix): 


F^{S-S) 


I E +YPn{n\S\n)\^, 

^ I -^mn _ 

m,n,m^n 


(16) 
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where for convenience the quantity 


^mn — Pm') — 2 ^n) 


is introduced. It can be shown by simple algebra that 

Fo{6S;6S)=Fo{S-S)-\{S)t\^ 

where 6S = S — {S)t- Since by definition 

{{6S‘^)'^)t := ^Pm(m|S'|m)^ - (S’)!, 

m 

where SS‘^ = S'^ — {S'^)t , we can present Eg. dTSl) in the form 


1 


i ^ \{m\S\n)\^^^^!^=Fo{5S-,SS)-{{SS<^)^)T. 


m,n,m^n 

In view of Eqs. m and (1201) from Eg. dldp . we obtain 

dlj,M = YFoiSS;SS). 

Recall that alternatively we have the well known relation (see, e.g. 0) 

'■i 1 d^\nZ{h) 


Fo{dS]SS)= [ dr {[dS{T)SS])j, = 
Jo 


/32 d^h 

Using the definition of the isothermal susceptibility with respect to the field h: 


|/i=0- 


Xh=o ■— 


1 In Z{h). 


the following result emerges 


P dK^ 


^2 _ /3,,t 

Up 


|/t=0, 


^BKM — -^Xh=0- 

We shall follow the above recipe in order to study the relations between different metrics. 


(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


B. Morozova-C'encov metric 


We shall show that the operator monotone function fMcix) = rUc ™^kes sense in physical applications. 

The corresponding Morozova-C'encov function 


CMciPm, Pn) = 


Pn In Pm \ Pn “i” Pm 
Pn Pm 


(25) 


yields the metric d\jQ already discussed in the Introduction. Since to our knowledge appropriate citation of this fact 
is unknown a simply deviation is presented. By inserting (1251) into m, we obtain 


— 7 ~~ + E \{''JT'\S\n)\‘^{pn + Pm)- 

^ Pm O 

m m,n,m^n 


(26) 


Now, using the relation 

{S^)T = \Y.\{m\S\n)\\pm + Pn) (27) 

m,n 

and Eqs. dH) and (IT^ we get the final result 

dlic = “ {F)t)'^)t- (28) 

Note that in the case of commuting operators T and S, according to Eo. d^ both metrics and dj^Q coincidece, 

as it must be. 
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C. Bures metric 

There are different ways to obtain the analytical expression of the Bures metric d\. It may be obtained from the 
infinitesimally close form of the Bures distance. For any two states pi and p 2 the Bures distance can be expressed in 
terms of the Uhhnann-Jozsa fidelity: 


T{pi,p 2 ) = Tr 



1/2 
P 2 P 1 , 


(29) 


as 


Db{pi,P2) = '\/2 - 2T{pi,p2). (30) 

Thus dg is obtained as the leading term in the expansion of Eg. (1301) . Among the operator monotone functions 
introduced by Morozova Cencov and Petz it is associated with the maximal operator monotone function fsix) = 

By setting it in Eq. m we obtain the result 


d 


2 _ 

B — 



dpi 

Pm 



m,n,m^n 


|(n|5'|m)p {pn-PmY 
|lnpm-lnp„p Pn + Pm 


(31) 


(see, e.g. and comments therein). 


D. Chernoff/Wigner-Yanase metric 


The quantum Chernoff metric is the infinitesimally close form of the Chernoff distance: 


Dqc{Pi, P2) — 1 — Q{pi,P2), 


(32) 


where 


Q{pi,P 2 ) = min Tr{pl V2) 

0<s<l 


(33) 


is known as the nonlogarithmic variety of the quantum Chernoff bound [ 2 ^ . It can be written as [1, [l^, [2l|, [ 2 ^ [s^ 



dpi 

Pm 



m,n,m^n 


\{ n \ S \ m )\^ 1/2 _ 1/2.2 

I 1 1 12 \Pm Pn / 

I Inprn - lnp„|^ 


(34) 


The metric (1341) coincides up to a factor 1/2 with the Wigner-Yanase metric m, i.e. (pQc = (l/2)d^y,, which may 

be obtained from Eq. m with the help of the operator monotone function fwvix) = . 

Here, the following comment is in order. In terms of Green’s functions a finite-temperature generalization of 
the fidelity susceptibility was proposed in [ 3 ^ . This quantity is quite different from the Bures-Uhlmann fidelity 
susceptibility (i.e. d^) and we denoted it as Xf- The only similarity between both metrics is that they have the same 
T = 0 limit. The fidelity susceptibility Xf is obtained as the first nonvanishing term in the expansion of the ’’fidelity” 


TT{pi.P2) = ^Tr{p\/^py^). (35) 

For for more details see ref. [l^. It is worth noting, that the spectral representation of y]) is presented by Eq. (14) in 
ref. [ 3 ^. It is simply to check that, in our notations, it emerges exactly as the quantum Chernoff metric Ea. (l34p . 
Since, the computational problems in studying y]) can be efficiently tackled by the quantum Monte Carlo approach 
this fact transfers new computational possibilities in the field concerning the quantum Chernoff metric. 


III. BKM AND MC METRICS AS REFERENCE METRICS 

In order to quantify the deviation of the one of monotone Riemannian metric d/ from the or we shall 

recast Eq. m in a form suitable for further elaborations. 
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A. BKM metric 


In order to discuss the deviation of a metric from the BKM metric, we need to rewrite Eq. (IT^ in the form: 

‘'? = jE^ + T E w (V-.) ~ f", K'”IS|")P, 


m,n,m^n 


In Pn - In Pr, 


(36) 


where we introduced the family of functions 


gf{x) := 


- 1 

2xf{e‘^^) ’ 


(37) 


These functions satisfy (due to the symmetry xf{x ^) = f{x)) the condition gf{x) = gf{—x) and g/(0) = 1. In the 
considered particular cases of /’s one is led to 


gHa 


, . sinh2a:: / n , / \ tanhia; tanhx , , 

(a:) = ,gBKM (X) = l,gwY (X) = — [X) = - ,gMC (x) = 


2x 


2^ 


tanh; 


(38) 


(We use argument X^n in the function g/(.) instead of ^ in the rhs of (1551) . remembering the relation Eg. dTTl) .') The 
presentation (1361) is convenient for two reasons. 

First, using the well known inequalities between hyperbolic functions one obtains 

gHar (X^n) ^ gMC (X^n) ^ PBKM (^mn) (— 1) ^ gWY (Xj^in) ^ PB (Xj^in) ■ (^^) 

From the above inequalities immediately follow corresponding inequalities between different metrics. 

Second, since functions g/ (Xmn) are even functions they have series expression with only even degrees of X^n- 
Thus, functions gf (X^n) can be expressed in a unified fashion as a formal series 


5/ (Xmn) — 1 + a 2 l-l{f)(.XjnnY^. 
1=1 


(40) 


We hope no confusion will arise by using the subscript 21 — 1 instead of 21 in order to enumerate the coefficients a’s 
in the above expansion. Particulary, in the case of Bures metric from the series expansion of the function x~^ tanhx 
it is easy to obtain 


a2i-i{B) 


( 2 ^ + 2 )! 


B 21 + 2 , 


(41) 


where B 21+2 are the Bernoulli numbers. The separation in coefficients ai{f) the very dependence on functions / allows 
to introduce the terms 


2 ^'a2i-i(/)(lnp„ - Inp^)^' - p„)|(m|S'|n)p (42) 

in the expansion Eg. (1361) . Remarkably, we shall show that these terms generate well known thermodynamic mean 
values. The presentation (1551) provides the role of Bogoliubo-Kubo-Mory metric as a reference metric. Recall that 
Eq. p6p simply presents the definition of choose gf (X^n) = 9bkm (Xmn) = 1- Thus gf (Xmn) 7^ 1 

measures the deviation of the corresponding metric dj from the Bogoliubo-Kubo-Mory metric. 


B. MC metric 


One may use an alternative presentation of instead of Eq. ([361): 

2 1 dp^ /3^ tv \ Pn + Pn 


■|(TO|S'|n)fy 


where 


m^n,m^n 


9f (x) ■■= 9f (a;) 


tanh a; 


(43) 


X 


(44) 
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Now, as a reference metric one can use (P^q. If we set / = Jmc in Eq.(j44]) then gMC i^mn) = 1 and Eq.([43|) reduces 
to the definition of Hence, functions gf serve as a measure of deviation of the corresponding d^’s from 

For the functions / presented in ([Mil , the application of definitions (l38l) in (l4^ yields 


, . 2(cosh2a:—1) ^ , tanhcc 

9Har[X) = - , -, gSKM (x) = - , 


(21)2 

ffMC (x) = 1, gWY(x) = 


2[1 —(cosha;) 


ffB(x) = 


tanhi 


n 2 


The even functions Ea. dTil) can be presented as a formal series 


(45) 


fff (^mn) — 1 + 0,2l{f){XmnY^ I 

1=1 


(46) 


which introduces in the summand of the expansion Eq.(j43]) the terms 

2 -(2i+i)a2;(/)(lnp„ - \np^Y\pn + Pm)\{m\S\n)\^. (47) 

These terms also can be expressed as some thermodynamical mean values (c.f. with Eg . (1421) '). In the next subsection 
we shall show that d\j(j appears as a first term of the obtained series representation of the monotone Riemannian 
metrics. 


C. Series presentations in terms of iterated commutators 


Let us consider the functionals 


ml 


(48) 


These functionals have been introduced in (see also the Appendix) as an useful tool to obtain different thermo¬ 
dynamic inequalities. By using the definition (H51) . one can present the terms ((1^ and (|47l) in the form 


and 


respectively. Now, the Ea. (l36l) reads 


Similarly, the Ea. (l4^ reads 


2-^^a2i-i{f)F2i{S;S) 

2-^^^+^^a2i{f)F2i+i{S; S), 

02 °° 

d) = 4iCM + ^ E 2-"'a2/-i(/)E2,(5; S). 

Z=1 


(49) 

(50) 

(51) 


«2 °° 

d)=dltcF^Yl 2-("'+')a2z(/)F2z+i(5; 5). (52) 

1^1 

Indeed, the first terms in Eqs. m and (|52ll are immediate consequences of the definitions (EID and (l28l) . while the 
summands follow from the terms and (IHH) . respectively. The key advantage of both formulas is that in the basis 
independent form functionals Fn{S; S) have the following presentations (see the Appendix) 

F„(5; S) = 2(-l)"+i/3"-i(i?„_iAo)T, (53) 

where the notion of iterated commutators 


Ro = Ro{S) = S, Ri = Ri{S) := [T, 5], ..., 

Rn=RniS) := [T,i?„_i(5)], n = 0,1,2,..., 


(54) 
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is introduced. (Note that nested commutators is also frequently used term.) It is worthwhile to emphasize the relation 
of Eq. (1551) with the moments of the dynamical structure factor Mp{S) (see the Appendix) 

E„(5;5) = 2/3’^-iM„_i(5), n = 0,1,2,... (55) 

With the help of Eg. (1551) one obtains from Eqs. (1511) and (1551) the result 

— ^BKM + i'o) (56) 

1=1 ^ ^ 

and 

rf/=rfMC + ^E(f) a2i{f)M2i{S). (57) 

1=1 ^ ^ 

This result attributes thermodynamical meaning to pure information theory ingredients. 

The series representations (1511) and (I52[) yield a proper definition of a monotone Riemannian metric provided the 
corresponding convergence condition is fulfilled. In order to remove the restrictions imposed by the convergence 
conditions one needs to perform an analytic extension of the obtained formulas. This issue in the particular case of 
the fidelity susceptibility (Bures metric) has been examined in the ref.fs^ by the examples of several popular models. 

The appearance of the iterated commutators (terms with n > 1) is a reminiscence of the well known Feynman’s 
disentangling procedure [13, . 

In the next Section we shall demonstrate that this allows the underlaying symmetry of the Hamiltonian to be 
efficiently explored and a closed form of the functionals (H51) to be obtained. 


IV. APPLICATION TO A MODEL 


Let us consider the Hamiltonian 1501 ^ : 


n{h) = kuj {qI -^+ hV^{Q+ + Q^), fc = 1,2,..., 


where Qt are operators obeying the commutation relations 


[Ql Q±] = ±Qf, [Qt, Qfc ] = MQl) - - 1), 


with the structure function 




(58) 


(59) 


(60) 


being a /c‘^-order polynomial in k. The Hamiltonian (1581) is employed in various physical problems (for definitions 
and a partial list of references, see (50l.[^b 

In this case the proposed approach is very effective since the iterative commutation between T = kuj {Q^ ~ ^) 

S = + Q^) implies some periodic operator structures after a finite number of steps [s^ 

i?„ = (-l)"R+=a”[Q^ + (-ira, ( 61 ) 


indicating an analytical expression as a function of n. The parameters k and ui enter in the c-number a = {koj)^^/l^ 
. Thus, the obtained series expansions Eqs. (1511) and (1521) can be used in a rather simple way to obtain closed-form 
expressions. 

The polynomial algebra of degree k — 1 defined by Eqs. dMl) has the following one-mode boson realization (50j |: 


Qt = 


Q-k = 




{■sfkf ’ {^fkf 

In terms of Eqs. (I62|) the Hamiltonian of the model takes the more familiar form 


'H{h)=ujh+h + h[{h+f+ w>0, fc = l,2,3,... 


(62) 


(63) 
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where bosonic operators 6, 6“*" obey the canonical commutation relations. 

It is worse noting that the k = 1 and k = 2 cases of (1551) give the Hamiltonians of the displaced and single-mode 
squeezed harmonic oscillators, respectively. The Hamiltonian (16311 for fc = 2 is also known as Lipkin-Meshkov-Glick 
(LMG) model in the Holstein-Primakoff single boson representation (see e.g. [l^ and refs, therein) and all the result 
obtained here can be related to this field. 

Inserting the expressions Ro and i? 2 rt-i in ([53|) . we obtain 

F 2 n{S;S) = -2ikPujf^-^IC{k), n = 0,1,2,..., k=l,2,..., (64) 

where 

JC(k)=k’^{[Q+-Q-][Q+ + Q-])T, k=l,2,... (65) 

Inserting expressions of Rq and i? 2 n in (1551) . we obtain 

F2n+i{S;S) = 2{kpujf^C{k), 

n = 0,l,2,..., k = l,2,..., (66) 

where 

£{k) = k'^{[Q+ + Q^]^)t, fc = l,2,... (67) 

Evaluation of the correlation functions Eqs. and dlZl) with the quadratic Hamiltonian T is now straightforward. 
The results for k = 1 and k = 2 are: 


/C(l) = -1, £{l)=2n + l, 

/C(2) =-2(2n-b 1), £{2)=4rL^, (68) 

where n = (e^"^ — 1)“^. With the help of Eas. dMl) and (I66p . Eas. (15T1) and (1551) may be recast in the form 


F — F 

“/ — “BKM 


/3^ f kpoj 


-1 r 


1-5/ 


/ kf^oj 


IC{k) 


and 


j2 _ ,2 

“/ — «MC ^ 


^ fk/3uj\ 


£ik). 


respectively. The relation (1551) on the particular example of Jb was found earlier in 
the whole class of the monotone Riemannian metric. 


(69) 

(70) 

Here we are able to consider 


V. SUMMARY AND CONCLUSIONS 


In this paper, we study model systems described by a Hamiltonian comprising the non commuting operators T 
and S via monotone Riemannian metrics. The formulas Eqs. dSB) and (l52l) are the key results of our study. They 
present the monotone Riemannian metric as a series in terms of the earlier introduced functionals Fn{S;S) (see 
also the Appendix). These are defined as a rather complicated double sum which however may be written in a 
basis independent form as thermodynamic mean values of n-times iterated commutators between T and S. The last 
provides significant computational advantage if the Hamiltonian is a linear form of the generators of some Lie algebra. 
In this case Fn{S; S) should be obtained in a closed-form expression as a function of n. Recall that the lowest-ordered 
functionals Fq{S;S) and Fi{S',S) are the Bogoliubov-Duhamel inner product and the symmetrized thermodynamic 
mean value of the operator S, respectively. Then Fq{6S;6S) is proportional to the BKM metric given by Eg. (1551) . 
while Fi{6S;6S) is proportional to the MC metric given by Eg. (1581) . These quantities are the starting point in our 
expansions Eqs. m and (1551) . 

From the linear response theory it is well known that the iterated commutators of an observable S with the 
Hamiltonian are related to the moments of the dynamic structure factor (DSF) through some sum rules This 

allows to present our results in an alternative form, Eqs. dMl) and dSZl). A look at the above formulas shows that 
the dependance on the operator monotone function / is only in the coefficients in front of the moments. If one aims 
to characterize a monotone metric on this setting, then moments of all orders (odd or even) should be considered 
together. 
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Formulas Eqs. (l56l) and m do not provide computational benefits per se. However they are very useful and 
informative because they transfer the information geometry problems into the realm of condensed mater physics with 
its wealth of methods for computing the DSF and its moments. For example, an instructive illustration may present the 
study where an application of appropriate (physical) approximations to Hamiltonian described one-component 
Coulomb plasma in thermodynamic equilibrium leads to explicit formulas for the arbitrary integer moments of the 
DSF expressed in terms of simple functions. 

The series representation (EH (or EH) bs regarded as a proper definition of the considered monotone Rie- 
mannian metric provided the corresponding convergence condition is fulfilled which can be checked on the framework 
of concrete models. 

We demonstrated our approach in an example with a Hamiltonian, expressed in terms of the generators of a 
polynomial deformation Lie algebra. Eg. (1581) . employed in various physical problems. It is shown that in this case, 
the infinite set of the moments of all orders (odd and even) can be written in a closed form, Eas. dMl) and Ea. dMl) . 
for the all monotone Riemannian metrics. Besides being of interest for its own sake the presented result may also be 
considered as a contribution to the linear response theory. 

At the end the following comment is in order. The implication of the immanent relations between the zero- 
temperature fidelity susceptibility and the moments of the DSF recently has been demonstrated in ref. [55j | . However 
the elucidation of the role of the moments Mp(S) (named as p-order generalized fidelity susceptibility) has been 
considered in quite different context. 
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APPENDIX: THE BOGOLIOBOV-DUHAMEL INNER PRODUGT, MOMENTS OF THE DYNAMIC 

STRUCTURE FACTOR AND SUM RULES 

Let us S is an arbitrary operator. The Bogoliubov- Duhamel inner product is defined us [29I - I3TI . [^ : 

{S;S)t= f'dXKssiX), (71) 

Jo 

where 

Kss{X)-.= {S+iX)S)^, S+{X)-.= e^P^S+e-^f^^, (72) 

and 


(•••)t := [Z(T)]-iTr{exp[-/3T]...} 

denotes the thermodynamic mean value. We warm the reader that definitions differ from EH by factor /3 and/or by 
involving on the first place in the Eg. (17^ the operator S instead of its adjoint S~^ exist in the literature. 

With the help of the Kubo identity (see, e.g.[^) 

[S,e-d'^] = -l3 [ dTe-d(^-^^'^[S,T]e-d^^, (73) 

Jo 

the following useful formula can be verified: 

I3{S;[T,S])t = {[S+,S])t, (74) 


where S and T are arbitrary operators. 

The definition of the Bogoliubov- Duhamel inner product [Eq. (1711) ] can be presented in terms of the spectral 
representation in which T is diagonal. Let us assume that the Hermitian operator T has a complete orthonormal set 
of eigenstates ]/) and eigenvalues Tf, T\l) —Ti\l). Then, in the basis of the eigenvectors of the Hamiltonian T we have 
the alternative form 


{S-S)t 





Pi Pm 
^ml 


\{m\S\l)\^ + Y,pms\l)\^ 

I 


(75) 
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where the notations {l\p\l) = pi and Xmi = 2 ^j3{Tm — Ti) are introduced. 

It is convenient to consider the spectral representation (for more details and history remarks, see e.g. 


/ OO 

dujQsiu 

-OO 


(76) 


where 


Qsiuj) = [Z{T)] ^^™|(n|5'|TO)p(5(a; - > 0, (77) 

m,n 

is the dynamic structure factor (DSF) relative to the operator S, where w is a real frequency and ujnm = Tn — Tm (h = 
1). DSF obeys the relation 

Qs(a;) = e^“Qs+(-w), (78) 


known as the principle of detailed balancing. 

Correspondingly the spectral representation of iFss(A) is given by Luttinger [si 

/ OO 

duiQsiuj; A), Qs(w; A) = Qsiuj)( 

-OO 

Hence, as follows from Eqs. (EU and (Hi we get 


,\pu: 


/ OO 

duj 

-OO 


-1 
(3uj 


Qs{uj), 


(79) 


(80) 


which is the original expression obtained by Bogoliubov 26j . For the further consideration it is useful to introduce the 
moments of the DSF 


/ + 00 

dioujPQsiio), P =-1,0,1,2,.... (81) 

-OO 

Then another alternative form of Eq. dn]) is 

(5;5)T = r'[M_i(5)+M_i(5+)]. (82) 

Eg. (15^ directly follows from Eqs. (1771) . (1751) and (15T|) . 

In our paper [s^ the functionals F„(S'; S'), n = 0,1, 2,... 

E„(S;S) :=2-i^|(m|S|l)npi-(-l)>„|.|X™ir-i, (83) 

ml 

have been introduced as a generalization of the Bogoliubov- Duhamel inner product [Ea. ((75l) ]. Some applications of 
these functionals have been demonstrated in |^, [^. Since Fb(S; S) = (S; S)t hereafter in the text we shall use 
this notation for the Bogoliubov-Duhamel inner product. Recall that Fi(S; S) = (S+S -f SS^)t- Using the notion of 
iterated commutators 

R^ = Ro{S)-.= S, Ri = Ri{S)-.= [T,S], ..., Ri = Ri{S) :=[T,Ri_^{S)]. (84) 

one has for even n = 21, Z = 0,1, 2, 3,..., 

F2i{S-S)= I3‘^\Ri-Ri)t = (85) 

where Rf denotes the Hermitian conjugate of Ri and by definition i?_i = Xst is a solution of the operator equation 

S=[T,Xst]. ( 86 ) 

In view of relations ((85)) . the functional F 2 i(S; S) may be called ’’Bogoliubov-Duhamel inner product of order I”. 

In the case of odd n = 21 + 1, Z = 0,1, 2, 3,..., one has 

F2i+i{S-S) = l3'^\[R+Ri+RiRi])T. 


(87) 
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It is remarkable that both formulas (|85l) and (I87|) may be written as one formula in the form 

+ n = 0,l,2,... (88) 

which is more useful in some cases. Eg . (155)) is obtained from Eg. dSSl) and Eg. (1571) with the successively using of the 
following identities 

{R^ Ri-i)t = {R^_^Ri)t, {Ri-iR^)t = {RiRI_i)t, 

{R^ Ri)t = {R^-iRi+i)t, {RiR^)t = {Rl+lRt-l)T, 

(89) 

which are simple consequence of the cyclic property of the trace and the definition of the iterated commutators (1841) . 
Setting n = 0 in Eq. (1551) one gets the relation 

Fo{S, S) = -r^[{R.iS+)T - {RtiS)T]. (90) 

It presents the Bogoliubov- Duhamel inner product as a thermodynamic mean value of the solution of Eg. (1861) and S. 

Erom the other side the functionals E’„(S'; S) defined by Eo. (l83l) may be obtained using the definition of the moments 
of the DSF Eo. ((5T]) and the relation ([75]), i.e. 

F„(5;^) =^"-i[M„_i(^) + M„_i(^+)], n = 0,1,2,... (91) 

Now, it is easy to check that 

M„_i(5) = (-l)"-i(i?„_i5+)r, M„_i(5+) = {Rt-iS)T, (92) 

which are not but the well known sum rules for the moments of the DSF in the linear response theory (see, e.g.fs^ 
and [s^, where some sum rules for the lowest moments are presented). The expressions Eg. (1571) provide an algebraic 
way to evaluate the moments of the DSF. 

The expressions (1551) and (1571) take a simpler form in the case S = S~^ that is required for the calculation in the 
text of the paper. From Eq. (1551) with the observation that = (—l)”i?„, n = 0,1, 2,... we obtain 

Fn{S;S) = 2{-ir+^f3^-\Rr,.iS)T, n = 0,1,2,... (93) 

or alternatively from Eq. m 


F„(5;^) = 2^"-iM„_i(,S), n = 0,l,2,... (94) 

We note that KssW has the properties of the scalar product and has been studied separately [H, [s^ because of 
its relation with the Wigner, Yanase and Dyson (WYD) skew information 

Is.sW = -l[Z{T)]-^Tr 5]) , 0 < A < 1. (95) 

via the relation 


KssW = Kss{0) - 7ss(A). 


(96) 


The Eq. (1551) allows to emphasize the relation of Fn{S; S) with the WYD information (1551) . In our further consideration 
we shall follow ref. [s^. Using the derivatives 


^ _^e-/3ATy^ ^ 

dX dX 

after some simple algebra one obtains 


(97) 


jn 

— iFss(A) = (-l)"/3"7^sfl„(A), n = 0,1,2,..., (98) 

where i?„ = [T, and i?o = S. The properties of these derivatives are studied in details in [sj emphasizing the 

relation with the Bogoliubov and Tyablicov Green’s function method. This seems to provide an useful method for 
their calculation. Our finding is that (if 5'+ = S) 


( 70 - 1 ) 

Fr.[S-S) = 2r'^^^if5s(A)U=o, 


n = 1,2,... 


(99) 
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or equivalently 


K(5;5) = -2/3-i^|j^/ss(A)U=o, n = 1,2,..., 


due to Ea. dMl) . Differentiating and integrating Eg. (17^ with respect to A, one obtains 


r 


—KssiX)=P^I duju:^Qsito-,X), n = 0,1,2,... 


and 


f dXKssiX) = /3 ^ f duju} ^Qsiuj;X), 

J J —OO 

respectively. At A = 0, these are proportional to the moments of the DSF 

/ OO 

duJOj'^Qsiuj), n =-1,0,1,2,... 

-OO 

Finally, let us note that from Eas. dMl) and (IIOII) we obtain Ea. dMl) as it must be. 
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